Elastic Convection in Vibrated Viscoplastic Fluids 
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We observe a new type of behavior in a shear thinning yield stress fluid: freestanding convec- 
tion rolls driven by vertical oscillation. The convection occurs without the constraint of container 
boundaries yet the diameter of the rolls is spontaneously selected for a wide range of parameters. 
The transition to the convecting state occurs without hysteresis when the amplitude of the plate 
acceleration exceeds a critical value. We find that a non-dimensional stress, the stress due to the 
inertia of the fluid normalized by the yield stress, governs the onset of the convective motion. 
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Complex fluids comprise a large number of fluids in 
everyday and industrial usage. Viscoelastic fluids often 
exhibit interesting behavior when subjected to mechan- 
ical stresses, such as rod climbing and the open siphon 
effect 0,0. Viscoplastic fluids that posses a yield stress 
can exhibit both elastic solid-like behavior and fluid-like 
behavior. Thus a blob of viscoplastic fluid can retain its 
shape when laid on a plate, even if turned upside down, 
if the stress due to its own weight does not exceed the 
yield stress. With the application of a mechanical shear 
stress larger than the yield stress, resistance to the stress 
abruptly decreases and the fluid starts to flow or falls in 
drops. In the present work, we observe a new instability 
that arises in soft yield stress fluids driven far from equi- 
librium that reveals characteristics particular to plastic 
flow which is not observed in various pattern formation 
in non- Newtonian fluids0,0| without yield stress. 

Vibrated materials often show a variety of patterns 
such as Faraday's surface waves in Newtonian fluids, 
stripe and labyrinthine patterns and convection of gran- 
ular materials^ 0, and localized excitation of shear- 
thickening fluids 7J. The convection of vibrated vis- 
coplastic fluids observed in the present work is a new 
example of pattern forming instabilities of vibrated ma- 
terials. In many pattern-forming systems, the evolution 
of the pattern can occur on a timescale much slower than 
that of the forcing. Understanding how the slow mode 
of the pattern dynamics emerges from the fast excitation 
mode is a fundamental problem in the nonlinear dynam- 
ics of far from equilibrium systems. 

When a viscoplastic fluid placed on a horizontal plate 
is vertically vibrated, freestanding convective rolls ap- 
pear as a lip shape, as shown in Fig. ??. We observe 
that a pair of rolls is spontaneously created, and counter- 
rotate with a slow circulation speed independent of the 
fast vibrating mode of the plate. The convection rolls 
are freely standing on the horizontal surface without side 
and top boundary walls. The diameter of the rolls is 
spontaneously selected and not determined by the depth 
of the container. This is in contrast to thermal convec- 




FIG. 1: (a) Overhead view of the convecting rolls on a hor- 
izontal plate (Rl gel), vertically vibrated at / = 40 Hz and 
acceleration amplitude V = 11.1. The volume of the gel used 
was V = 6 mL. (b) Resulting roll pattern using a smaller vol- 
ume, V = 3 mL, at / = 40 Hz, F = 12.0. The scale bar is 
the same as that of (a), (c) Ring pattern obtained at higher 
volumes. The mean convecting flow occurs in the outer ring 
region. / = 40 Hz, F = 15.0, and V = 9 mL. (d) Rolls with a 
branch can be formed with a different initial condition. Os- 
cillation parameters are / = 100 Hz, V — 32.6, and V — 5 
mL. Arrows represent the direction of the flow. 



tion where the depth of the fluid/container constrains the 
wavelength of the convection. The critical condition for 
the onset of convection is independent of the volume and 
the height of the fluid for a wide range of parameters, 
and is determined by a newly denned parameter: a non- 
dimensional stress produced by the inertia of the fluid 
interacting with the oscillating plate. 

The experimental apparatus consists of a horizontal 
aluminum plate (120 mm diameter), mounted on a vi- 
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bration device (VG-100C: Vibration Test System). The 
plate is vertically vibrated with a sinusoidal oscillation 
with amplitude A and frequency /, i.e. z — Asin(27r/t), 
by a function generator (NF 3540). The dimensionless 
acceleration amplitude is defined as T = A(2irf) 2 /g, 
where g is the gravitational acceleration. The accelera- 
tion amplitude of the plate is measured by an accelerom- 
eter (PV-90B:RION) to an accuracy of ±0.2 m/s 2 . The 
motion is recorded by a high speed CMOS camera with 
a 1024 x 1024 pixel resolution and time resolution of 
l/1000s. We used fluids with a monotonically shear- 
thinning property above the yield stress strength, i.e. 
viscoplastic. We report the results for a typical shear- 
thinning fluid gel (sol-gel intermediate), commonly used 
for ultrasound examinations (Dane-gel Rl, Rohde prod- 
ucts). The material is made of a polymeric gelling agent 
and water. However, we observed qualitatively similar re- 
sults using other shear thinning yield stress fluids such as 
xanthan gum, shaving gel, and toothpaste. We measured 
the viscosity of the material using a cone plate rheome- 
ter with temperature control (Rheosol G5000, UBM). 
We confirmed that the viscosity decreases with increas- 
ing shear rate. We estimated the yield stress of the gel 
as tb — 150 Pa, using the same device by measuring 
the maximum stress supported by the gel under constant 
shear. 

A fixed volume of the gel is placed on the plate. The 
initial shape of the gel is amorphous. The gel is elas- 
tically stretched and compressed during the plate oscil- 
lation period. As the plate acceleration is increased it 
becomes apparent that the gel is fluidized, and it rapidly 
collapses and spreads onto the plate. A lump initially 
laid on the plate is thus gradually encircled by a rim of 
inward rotating rolls as the acceleration is increased. The 
rim of roll elongates in one direction along the horizontal 
surface breaking the azimuthal symmetry. Within a few 
seconds the middle part of the lump is invaded by rolls, 




FIG. 2: Side view of the oscillating motion of the rolls at 
various phases of the plate oscillation (/ = 40 Hz, F = 14.8, 
and V = 5 mL). Time proceeds from left to right at 5 ms 
intervals. The oscillation of the gel is locked to the plate 
frequency but out of phase [see Fig. |S[c)]. However, the 
mean convective flow occurs on a time scale independent of 
that of the oscillation period. 
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FIG. 3: State diagram of the materials as a function of accel- 
eration and frequency. The fitted line shows the state bound- 
ary. Below the line the rolls do not rotate. Above the line 
they can rotate if the roll shape has been clearly formed. The 
dashed line is obtained by a linear fit. Inset: The same plot 
for non-dimensional stress defined by Eq. Q. 



and finally become lip-shaped as shown in Fig. ^(a-b). 

Detailed dynamics can be seen with the use of a high 
speed camera. The rolls oscillate elastically at the same 
frequency as the plate (shown in Fig. [2J) but with a dif- 
ferent phase [shown in Fig. [3(a)]. The gel appears to be 
folding with each oscillation of the plate. Each roll ex- 
hibits inward rotation: outer edges lift from the surface 
of the plate and fold toward the center of the lump to 
create rolls as shown in Fig. [21 This motion is very pe- 
riodic and therefore the net displacement of each part of 
the fluid is small after each complete period of the oscil- 
lation. However, we observe a relatively slow circulating 
convection within the rolls, by tracking the motions of 
tracers synchronously with the period of the vertical os- 
cillation. The speed of the convective motion (e.g. ~ 1 
Hz) is much slower than the fast mode of the periodic de- 
formation due to the plate oscillation (e.g. 40 Hz). The 
motion is independent of the frequency of oscillation but 
increases as a function of the acceleration of the plate 
above the onset. 

We measure the onset acceleration amplitude r c to de- 
termine the dynamical state diagram plotted as a func- 
tion of frequency, shown in Fig. Once the roll pattern 
is formed, the shape of the rolls persists even after the 
vertical vibration is stopped due to the high viscosity and 
yield stress of the gel. If we decrease T from a convecting 
roll state, rotation of the rolls ceases immediately when 
r becomes lower than the critical value r c and starts as 
we increase T above T c without hysteresis for the onset. 
The state boundary falls on a line i.e., T ~ (0.27±0.01)/. 
Therefore the convection starts when the velocity ampli- 



3 



■a 

2 6 



as 

£ 2 



? 6 
84 

I 2 

< 



5 10 15 20 

Acceleration r 



r' 1 Y A ' 

-J ^r- 

I'f JiT* --40Hz 

5-*-'* ^-60Hz 



2 3 4 

Acceleration r-T 



FIG. 4: The mean angular velocity Q as a function of the 

acceleration above onset, F — r c . Results for / = 40 Hz with 

V = 5.0 mL and 60 Hz with V — 2.2 mL. The horizontal axis 

is shifted such that the onset value r c is at the origin. Each 

is fit by f2 = a(F — Fc) 13 , where a and /3 are fit parameters. 

This results in Q = (0.82 ± 0.04) X (T - Pc) 1 - 50 * - 04 [s _1 1 and 

M. 55 ±o.03 r s -ii Inset 



n = (0.47 ± 0.02) x (r - r c 

plot without subtracting r c . 
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tude v = A(2tt/) is larger than a certain threshold ve- 
locity v c = 0.39m/s. This threshold varies with different 
kinds of testing materials. Therefore it is reasonable to 
introduce non-dimensional parameters in order to clar- 
ify the mechanism of the onset of convection. Here we 
introduce a non-dimensional parameter, a stress exerted 
by the plate due to the inertia of the material normalized 
by the yield stress cr v , 



M 2 (2^/) 2 



(1) 



The stress due to inertia is estimated to be a = 
pA ■ A(2irf) 2 — pv 2 , where p is the mass density of the 
material, £ is the ratio of the exerted stress and the yield 
stress. We plot the non-dimensional stress £ for vari- 
ous volumes of material and various frequencies, in the 
inset of Fig. The critical condition for the onset of 
convection is £ ~ 0(1). Despite the fact that the esti- 
mation of the yield stress has some ambiguity and the 
stress distribution is not uniform, the data for various 
different conditions cluster near the line corresponding 
to £ = 0.8. The result strongly supports the hypothesis 
that the convection starts when the stress distribution in 
the material exceeds the yield stress a > <j y , i.e. above 
the solid to fluid transition. 

We next measure how the mean angular velocity of the 
convection rolls changes above the onset. Rotation of the 
rolls is monitored by following small bubbles or density 
matched tracer particles which are immersed in the gel 
for visualization. The convection speed of the fluid de- 
pends on the location within the convecting rolls. The 



rolls rotate faster in the middle part and slower at the 
ends. The mean angular velocity (radian/sec) near the 
mid point of the rolls is obtained from the measurement 
of the time interval required for several complete cycles of 
the tracers. The result is shown in Fig. The errorbars 
in Fig. 0] are relatively large due to the intrinsic angu- 
lar velocity distributions in the rolls as well as statistical 
variation. We find that the mean angular velocity in- 
creases as an increasing function of the plate acceleration 
T above the onset and is proportional to (r-Tc) 1 ' 53 * 03 . 

Now we note an important feature of the convecting 
rolls: The diameter, not the length, of the convective 
rolls is spontaneously selected. This is particularly true 
when the frequency of the plate is relatively low. i.e. 
30Hz - 40Hz. Figure Ufa) shows the time evolution of 
the diameter of one roll observed from above using the 
high speed camera. We obtain Fig. (b-c) by plotting 
the maximum and the minimum observed roll diameter 
dependence on acceleration and volume. The maximum 
and the minimum values for the diameter of the rolls are 
independent of acceleration and volume well above the 
onset, and for sufficient volume of the material at lower 
frequencies, e.g. / ~ 40 Hz. For higher frequencies, 
e.g. / ~ 80 Hz, the diameter of the rolls is affected by 
factors such as the initial conditions, suggesting that the 
stability band of the wavelength selection is wider for 
higher frequency. We can observe very narrow rolls for 
the same volume of the material compared with that at 
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FIG. 5: (a) Dotted line: Time evolution of the diameter of 
a single roll as observed from above, looking down onto the 
plate, for / = 40 Hz, V = 5 mL, and F = 14.3. Solid line: 
Position of the plate showing the phase difference, (b) De- 
pendence of the range of roll diameters on acceleration V for 
/ = 40 Hz, V = 5 mL (F c = 10.1). The upper (lower) dots 
show the maximum (minimum) value of the measured diam- 
eter. Note that there is a plateau region for the maximum of 
the diameters when the acceleration F is well above the onset 
value Fc- (c) Size dependence of the roll on the volume of 
the material measured for / = 40 Hz, F = 12.3 (above Fc). 
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FIG. 6: Aggregation of the gel in a vertically vibrated cylin- 
drical container (acrylic resin walls and aluminum plates). 
Time proceeds from left to right at intervals of 2 s. The gel 
eventually gathers into one area enclosed by rolls and the 
wall, vacating the other areas. The experimental parameters 
are fixed at / = 50 Hz, V = 27.3, and V = 25 mL. 

low frequencies. We also note that pattern selection has 
some initial condition dependence as is typical of other 
pattern forming phenomena far from equilibrium. For 
example, a triangular vertex of roll pair or a concentric 
ring of roll were observed as shown in Fig. ^c-d), by 
imposing different initial conditions. 

The measurements described until this point were con- 
ducted on a freestanding blob of gel on a horizontal plate 
without sidewalls. However, convection of the gel still 
occurs in the presence of sidewalls. When placed in a 
container, the gel self-segregates into a "material sea" 
region and a "dry" region, with convection rolls forming 
the "shores". This process is shown in Fig. H3 An ini- 
tially uniform shallow (depth of about 15-20 mm) layer 
of the gel was placed covering the bottom plate of the 
cylinder. As the plate acceleration is increased above 
the critical acceleration, a roll eventually forms near the 
wall. Once a roll forms, it aligns itself perpendicular to 
the point of contact with the wall (see Fig. . The roll 
tends to straighten and thus the action of the roll is to 
peel the gel from the bottom plate and move it away from 
the wall. A very thin layer of gel remains stuck to the 
plate in the vacated region. Eventually the gel is thus 
segregated into two regions: a shallow region bounded 
by the container wall and rolls, and a much thinner layer 
of gel in the region vacated by the rolls. 

The loss of stability and subsequent transition to the 
convecting state of a vertically oscillated viscoplastic fluid 
occur due to the solid to fluid transition, as evidenced by 
the order unity transition of the nondimensional stress. 
Above the onset, we suspect that a nonlinear mecha- 
nism such as acoustic streaming may be responsible for 
the emergence of the slow mode from the fast oscillat- 
ing modeja, It is reported that the circulating mo- 
tion around a vibrating rod in viscoelastic materials is 
in the opposite direction compared to that in Newto- 



nian fluids. The fact that the flow on the center line 
is directed toward the rod coincides with our observa- 
tion. However, the observed velocity of our convection 
and that of acoustic streaming behave very differently as 
a function of vibration velocity. We also speculate that 
the anisotropy of viscosity plays a crucial role in this 
phenomenon. Decrease of viscosity (shear thinning) de- 
pends on the shear stress tensor. Since the distribution 
of the stress is anisotropic, the viscosity coefficients of 
the fluid are also anisotropic. Once the convection rolls 
are formed by symmetry breaking, the shear in the az- 
imuthal direction is larger than the longitudinal direc- 
tion. This results in a smaller viscosity in azimuthal 
direction compared with the longitudinal one, and this 
effect may stabilize the motions along the azimuthal di- 
rections and straighten the rolls. 

The newly found state, freestanding convection, can 
become a new demonstration for viscoplasticity, the pres- 
ence of a yield stress in a complex fluid. Similar phenom- 
ena can be expected to occur in materials ranging from 
yield stress fluids to amorphous solids pH lill hj and 
glassy materials|l3t fl4| . Furthermore, enhanced mixing 
due to stretching and folding of the fluids during convec- 
tion may allow for applications to mixing of such fluids. 
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